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ABSTRACT 

Gamma-ray bursters emit a small fraction of their flux in X rays, and because X-ray 
detectors are often very sensitive they may probe the gamma-ray burst universe more 
deeply than the current best gamma-ray instruments. On the reasonable assumptions 
that spectra of bursts observed by BATSE may be used to predict the X-ray fluxes 
of gamma-ray bursts, and that any corona of bursts around M31 is similar to the one 
around the Milky Way, we predict the rate at which the wide field cameras on board 
BeppoSAX should detect bursts from the Milky Way and M31. These rates are such 
that a one-month observation of M31 would have to either detect bursts from M31 or 
exclude most galactic models of gamma-ray bursts. (It is shown how the remainder can 
be dealt with.) Therefore such an observation would settle the long-standing dispute 
over their location. 
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1 INTRODUCTION 

The results of the BATSE mission (see Fishman and Meegan 
1995), combined with earlier data sets for bright bursts such 
as the one collected by PVO (Fenimore et al. 1993) have 
shown that (1) gamma-ray burst positions are distributed 
uniformly and randomly on the sky (Briggs et al. 1996a,b) 
and (2) the cumulative number as a function of peak flux, 
A(> By), is consistent with a constant rate density of bursts 
within some volume around us, and a decreasing density out¬ 
side that volume. This implies that we are at the centre of 
a gamma-ray burst universe of which we can see the edge 
and which looks the same in all directions. Most distance 
scales are therefore excluded. The first remaining one is the 
high-redshift universe, with the edge being caused either by 
cosmological volume effects near and beyond z = 1 or by 
evolution of the density at moderate redshift (or both). The 
second one is an extended corona of our Galaxy, much big¬ 
ger than the dark-matter halo and invented for the purpose 
of housing gamma-ray bursts (GRB). We are not strictly in 
its centre, but the average GRB distance can be made large 
enough that the anisotropy due to our offset from the centre 
is below the limit set by the BATSE data on burst positions. 
At the same time it can still be small enough that we need 
not see M31 (Briggs et al. 1996b). The aim of this paper 
is to demonstrate the capability of the Wide Field Gameras 
(WFC) on board BeppoSAX (launched in April 1996) to 


distinguish between these options by searching for the hy¬ 
pothetical corona of GRB around M31. We first discuss the 
Z-ray detectability of GRB (Sect. ^) and our implementa¬ 
tion of corona models (Sect. H). Our results are presented in 
Sect. ^ and compared to previous results in Sect. 

2 X-RAY DETECTION OF GRB 

Ginga observations and some earlier detections indicate that 
gamma-ray bursts emit some X rays (for an overview of early 
X-ray detections, see Preece et al 1996). It is only a small 
fraction of the flux (2% or so median; Laros et al. 1984, 
Yoshida et al. 1989), but since X- and gamma-ray instru¬ 
ments are photon counting devices it is the higher count 
rate ratio of X rays to gamma rays that matters. In addition. 
X-ray detectors usually have lower backgrounds because of 
their imaging capability. In all, they can see some fainter 
bursts than BATSE, the currently most sensitive gamma- 
ray detector looking for GRB, at a price of having a much 
smaller field of view. This was recently used by Hamilton 
et al. (1996) to constrain galactic-coronal models of GRB us¬ 
ing archival Einstein data. They used galaxies typically a few 
Mpc away, the compromise distance for Einstein’s exquisite 
sensitivity but very small field of view. The lower sensitivity 
of the WFC mean we should observe more nearby galaxies 
to constrain galactic-coronal models with them; the large 
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field of view means that we can go as close as M31 without 
losing too many bursts because they lie outside the field of 
view. 

2.1 The WFCs on board SAX 

There are two WFCs on board SAX, looking at opposite 
directions on the sky and perpendicular to the on-axis in¬ 
struments. They are coded-mask imaging instruments with 
an entry mask of 256 x 256 1 mm^ pixels placed 700 mm away 
from the detector plane. This leads to a response that is ap¬ 
proximately triangular in both x and y and falls to zero 
20° away from the optical axis along the x and y direc¬ 
tions. It is sensitive to the energy range 1.6-32 keV. The re¬ 
sponse function was kindly supplied to us by Dr. J. Heise of 
SRON Utrecht. The angular resolution is a few arcminutes 
for bright sources. The 130 cts/s background of the instru¬ 
ment is mainly due to the diffuse extragalactic emission inte¬ 
grated over the field of view. The instrumental backgrounds 
are small and stable due to the low equatorial orbit which 
avoids the radiation belts and the South Atlantic Anomaly. 

2.2 The X-ray fluxes of GRB 

Previous workers have used a mean flux ratio of typically 2% 
between the BATSE flux and the X-ray flux of a gamma-ray 
burst. Rather than rely on the few X-ray detected gamma- 
ray bursts, we note that the WFC and BATSE sensitiv¬ 
ity ranges overlap in the 10-30 keV range, and that there¬ 
fore extrapolating the BATSE spectra into the WFC band 
should give reasonably good estimates of the expected SAX 
WFC count rates. Band et al. (1993) published detailed 
spectral fits to a sample of bright GRB from the first-year 
BATSE catalogue. The model consists of two power laws 
connected by a smooth transition at a break energy. Almost 
all break energies are well within the BATSE range. The 
model fits are shown in Fig. scaled to a photon number 
flux of 1 phot cm“^ s”*^ integrated over the 50-300 keV band. 
(The integrated count rate over this band is used by BATSE 
to trigger bursts.) We folded each of their 54 best-fit models 
through the WFC response matrix and computed the ra¬ 
tio $wFC of the WFC on-axis count rate to the 50-300 keV 
photon number flux. This may seem an odd ratio to take, 
but the WFC detectability is determined by the count rate, 
whereas the GRB flux distribution seen by BATSE is usu¬ 
ally reported after correction for detector response, i.e. as 
a photon number flux. A histogram of these ratios, and of 
the more commonly used F{2 — 10keV)/F'(50 — 300keV), is 
shown in Fig. The median flux ratio agrees with the Ginga 
estimates of 2-6% (Yoshida et al. 1989). But the large spread 
is crucial because it implies that a substantial fraction of 
GRB will be brighter than previous authors have estimated 
using a constant Fx/F-y and will therefore be more easily 
detected. 

A possible concern is the extent to which our sample 
is representative of the population. The selection criterion 
employed by Band et al. (1993) to obtain their subsample of 
the first BATSE catalogue for spectral fitting is essentially 
a brightness cut, and the authors consider the sample to 
be effectively a complete flux-limited one. This reduces the 
question of fair sampling to whether the spectral shape sys¬ 
tematically changes with flux. There is one known effect of 



E [keV] 

Figure 1. The best-fit model spectra for the GRB sample from 
Band et al. with their extrapolation into the WFC band. Note 
that the range over which the model fits were made significantly 
overlaps with the WFC band. 

this kind, namely that faint bursts have break energies that 
are smaller by a factor 2 than bright ones (Malozzi et al. 
1995). Since our spectral sample is from bright bursts, we 
may calculate the effect on the values of 'hwFC by decreas¬ 
ing the value of the break flux by a factor 2 for the whole 
sample. Each value of 'I'wfc changes differently because the 
change depends on the spectral slopes, but the resulting dis¬ 
tribution of 'hwFC is shifted to higher values by a factor 2, 
i.e. faint bursts are relatively brighter in X rays. 

Moreover, if there is a correlation between spectral 
shape and luminosity, and this is translated into a weak de¬ 
pendence of spectral shape on flux via the selection bias on 
flux, we probably should not even want to correct it because 
we will only see the most luminous bursts from M31 with 
SAX, so a sample of spectra that favours luminous bursts is 
a better one to use. Similarly, since the sample is gamma- 
ray selected, any bias in the spectral sample will be towards 
GRB that are gamma-ray bright, so the true population is 
likely to have on average greater X-ray brightnesses if the 
sample is not fair. 

In what follows we shall use the 'hwFC distribution from 
the Band et al. sample to estimate the X-ray detectability 
of GRB. We note that all the known and potential biases 
discussed would increase the number of detectable bursts 
from M31 over the calculations presented below. 

2.3 X-ray detectability given a gamma-ray flux 

To decide whether a given GRB will be detected by the 
WFC, we assume that it is a standard candle in the BATSE 
band. It has already been shown (Hakkila et al. 1995) that 
the range of gamma-ray luminosities of GRB must be small 
for all corona models that are still viable. Ulmer and Wijers 
(1995) also showed that for most luminosity functions, the 
luminosity distribution of detected GRB is narrow even if 
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Figure 2. The distribution of classical X-ray flux (2-10 keV) to 
BATSE (50-300 keV) flux ratio (top) and the WFC (1.6-32 keV) 
to BATSE photon number flux ratio, which is more relevant for 
burst detection. The large difference is partly due to the energy 
per photon being much less in X rays, and partly due to the fact 
that the WFC band is much wider than the classical X-ray band. 


the intrinsic luminosity function of the population is not. 
If the luminosity function is wide it increases the spread in 
X-ray luminosities at constant mean and therefore helps X- 
ray detection, unless there is a strong correlation between 
luminosity and spectral slope, so the standard-candle as¬ 
sumption will yield the lowest estimate of the number of 
X-ray detectable GRB. 

The procedure to decide whether a model GRB will 
be detected by the WFC is then as follows. First, the 50- 
300 keV peak flux is given by the model. Then the WFC 
on-axis count rate is computed by multiplying with $wfc. 
The model also specifies the location of the burst, from 
which we calculate the position (x, y) of its image in the 
detector plane (to be precise, the intersection of the line 
connecting the burst location and the centre of the mask 
with the detector plane). The count rate due to a burst at 
(*, y) is less than that due to one on axis by a factor R{x, y), 
which accounts for the fact that it only illuminates a frac¬ 
tion of the detector and for the usual factor cos 6 to account 
for the fact that the detector plane is not perpendicular to 
the direction to the burst: 

= , 1 ) 

{x and y are measured in mm from an origin at the detec¬ 
tor centre, and the x and y axes are parallel to the edges 
of the square mask.) Given an integration time T and back¬ 
ground count rate b, we get the total number of source counts 
S — fP~f^wFcRT and background counts B = bT. The fac¬ 
tor / is required because P-, is the peak flux, which will gen¬ 
erally not be sustained for the full time T. In terms of the 


instantaneous number flux p^{t) from the burst (assumed 0 
outside the interval (0, T)) we can formally define / as 


/ 


Jo P-iiD df 

P^T 


( 2 ) 


Obviously, for a fixed value of T we will find a different / for 
each burst and / will usually be smaller for shorter bursts. 
We assume the background noise is Poissonian and large 
enough to be approximated by a normal distribution, so we 
can express the requirement that the burst be more than a 
standard deviations above the background as a constraint 
on P^: 


P-y > Pj 


ay/b 

f^wFcR{x, y)'/T 


(3) 


So we have now phrased the X-ray detectability as a con¬ 
straint on the gamma-ray flux, which is convenient because 
most of the modelling of GRB populations is done in terms 
of the latter quantity. To account for the fact that i&wFC has 
a distribution of values rather than a fixed value, we treat a 
burst at a given {x, y) with a given P~f as being detectable 
with probability 


i=l 


(4) 


where Six) is the Heaviside step function and we have de¬ 
fined P^i’^ as the minimum detectable flux for the ith sam¬ 
ple member. In other words, we just add up the fraction of 
model spectral shapes in the Band set for which its X-ray 
flux is above threshold. Since we do Monte Carlo simulations 
to find the rate of GRB detection (Sect.^ we can then sim¬ 
ply add up the Pdet values in each sky and flux bin to get 
the detected rates. 

In principle, the signal-to-noise of an off-centre source 
will be higher once it has been located, because the back¬ 
ground need only be taken over the part of the detector that 
is illuminated by the source. However, we will not be aware 
of such sources until they are first noted in the full data 
stream so this does not change the detection rate. 

For the preliminary investigation in this paper, we shall 
use T = 20 s and / = 0.5. This means we concentrate on 
the long part of the bimodal duration distribution of the 
bursts, which contains about 80% of the ones detected by 
BATSE. More advanced search techniques are clearly pos¬ 
sible. For example, one can use a number of trial values of 
T or even trial sky positions for the bursts to enhance the 
sensitivity of the search. At the same time, this would in¬ 
crease the number of attempts at detection and therefore 
require a higher signal-to-noise threshold to avoid spurious 
signals. It is not meaningful in our view to explore these 
possibilities here because the optimal strategy will depend 
on details of the data set we will obtain. There will be real 
signals from other sources, such as X-ray bursts and flare 
stars, that may well constitute a higher contaminating rate 
than the expected foreground of Galactic GRB. A good frac¬ 
tion of those should be discernible from GRB on spectral or 
temporal grounds, but near the threshold complete subtrac¬ 
tion may not be possible. Also, there could be instrumen¬ 
tal effects that very occasionally give a spurious signal, and 
the importance thereof and consequences for the search and 
analysis strategy are hard to predict. 
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In Fig. 1^ we show the detection probability for bursts at 
the centre of the detector field of view as a function of P^. 
The thick dashed curve is for BATSE, and the thick solid 
one is for the WFC using the nominal Band spectral set. 
Comparison of the two shows that above 0.3 phot cm~^ s“^ 
BATSE is more sensitive, but the curve for the WFC extends 
to well below the absolute limit for BATSE. This advantage 
vanishes if we look at bursts far away from the detector 
centre: 10° away along the detector diagonal the curve shifts 
up in flux by a factor 2.4, and 20° away that factor has grown 
to 12. So the effective field of view over which the WFC can 
probe fainter bursts than BATSE has a radius of about 10°. 

There is one more spectral effect, however, that could 
dramatically increase the detected rate by the WFC: in a 
recent paper, Preece et al. (1996) find that many GRB have 
an excess flux above the Band et al. fits at energies between 
5 and lOkeV. Inspecting their figure 5, it appears that about 
half of all GRB have an excess that is a factor 10 at 5keV 
and decreases to near zero above 10 keV. To investigate the 
possible effect of this excess, we model it as a multiplicative 
factor Cpreece that is constant below 5keV and decreases 
smoothly from 10 to 1 between 5 and 12keV. To be precise, 
we multiplied the spectrum by 



Cpieece{E) = 10 


= 1 + 4.5(1 + cos 

= 1 


E <5 
■k{E - 5) , 
7 

12 < f; 


5 < 


(with E in keV.) Adding this excess to each of the 54 spec¬ 
tra and recalculating the detection probabilities, we find that 
the entire sensitivity curve shifts down in flux by a factor 
of 7 (thin dashed curve in Fig.^. Since only about half the 
bursts may have this excess, a more realistic case would be 
an even mixture of bursts with and without excess (thin 
solid curve). We stress that the fraction, /ex, of bursts with 
excesses and the form of the excess are quite uncertain, so 
we explore the range 0.1 < /ex <0.9 here. But since we will 
have spectra of all detected bursts in the eventual observa¬ 
tions we can consistently account for it in the data, since 
even if no excess of GRB from M31 is found, a population 
of bursts with X-ray excesses would also greatly increase the 
foreground rate of GRB from our own Galaxy. Therefore we 
are in no danger of eventually overestimating the constraints 
on halo models from WFC data. 


3 CORONA MODELS 

For standard-candle gamma-ray bursts, the rate density of 
observable bursts as a function of distance from the centre 
of the corona follows directly from the observed N{> P). It 
is usually approximated as 


p{R) = 


po 


1 + (R/Rc)“ 


( 6 ) 


and has three parameters, the central rate density po, the 
core radius Rc, and the exponent a. Since the best values 
of a are small enough that the integrated density does not 
converge at large R, we have to add as a fourth parameter 
a cutoff radius of the corona, Ro, beyond which p = 0. An 
additional model parameter is the standard-candle gamma- 
ray photon emission rate . It turns out that we only need 


Figure 3. The detection probability of GRB at the centre of the 
WFC field of view compared with the BATSE trigger efficiency 
(thick dashed curve). The thick solid curve is for the Band et al. 
models; the thin dashed curve is with each burst given an X-ray 
E <ek?ess suggested by the data from Preece et al. (see text). The 
thin solid curve is for a mixed population of 50% bursts with an 
X-ray excess and 50% without. The integration time for the WFC 
curves is 20s, the background rate 130ctss“^, and the detection 
threshold is set at 5(t. 


two of these, Rc and Ro, because the others follow from 
them if we use known observational constraints. Moreover, 
the results depend only weakly on Ro. We briefly indicate 
how a complete model is defined once the two radii are given: 
First, we note that the break in the counts slope from —1.5 to 
a smaller value occurs at Rbreak = 20 phot cm“^ s“^. Since a 
burst with this flux is at distance Rc, this fixes the standard- 
candle value as = drrRoRbreak. Next we use the fact that 
BATSE observes 300 bursts per year per 47r steradians above 
Rcomp = 1 phot cm“^ s“^. Because this flux is 20 times less 
than the break flux, it corresponds to a distance RcVW, and 
we fix Po by requiring the integrated rate up to that distance 
to be 300/yr. a follows from the fact that N{> PA oc 
at the faint end of the BATSE distribution. Because the 
asymptotic counts slope at low fluxes for Eq. (^ is (3 —a)/2, 
this implies a = 1.6. This completes the model for given Rc 
and Ro. 

The reasonable range of Rc and Ro to explore is also 
limited by data. First, for Rc ;$ 30 — 40kpc BATSE would 
have detected an anisotropy due to our offset from the Galac¬ 
tic centre (Briggs et al. 1996b), and for Rc ^ 70kpc it would 
have seen M31. A minimum value for the outer radius follows 
from the fact that BATSE sees no sign of a truncation down 
to its 50% completeness limit of 0.28 phot cm“^ s”*^. The dis¬ 
tance to bursts of this flux is 8.5Rc, i.e. at least 250 kpc, so 
we conclude that Ro ^ 250 kpc. Due to the low sensitivity 
of our results to Ro we fix it at half the distance to M31 in 
our calculations, noting that this is by no means required by 
all models. (For neutron-star ejection models with beaming, 
we use Ro = 2Mpc; see below.) 
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Figure 4. Four images of the 
rate of detectable bursts from 
M31. Contours at at 10, 30, 
50, 70, and 90% of the peak 
rate, which is indicated as ‘max’ 
in each panel (in units of 
deg“^yr“^). The x and y are 
offsets from the detector centre 
in degrees. All maps are made 
assuming 50% of GRB have an 
X-ray excess. 


As stated, all standard-candle galactic-corona models 
have to satisfy this model for the effective rate density. Since 
non-standard-candle models have greater difficulty satisfy¬ 
ing the BATSE constraints and also lead to easier detection 
of M31 we shall not consider them here. However, the net 
rate density is obtained in very different ways by different 
models, and this has some effect on detecting M31. It does 
not matter whether a static halo is used or one in which the 
bursters fly out of the Galaxy at high velocity, because all 
models of the latter type contain a provision of gradual or 
delayed turn-on of the bursting mechanism to ensure that 
the net rate density becomes the same again. However, some 
models (e.g. Duncan et al. 1993) invoke beaming of emission 
from fast neutron stars along their velocity direction to avoid 
anisotropy. The opening angle of the beams required in such 
models is of order the distance Rgc from us to the Galactic 
centre divided by Re- This does matter, because it means 
that only bursters in M31 that move roughly towards or 
away from us will be seen. While this reduces the number of 
visible bursts, it also limits region of the sky where they are 
seen to a circle of angular radius about R^smO/D around 
the centre of M31, where 6 is the opening angle of the beam¬ 
ing cone and D is the distance to M31. Since 9 ~ Rgc/Rc 
the angular radius is Rgc/D, independent of core radius. 
This area is quite small (typically a few square degrees) and 
therefore the background in it is very small, which may com¬ 
pensate for the lower expected rate to still give a detectable 
excess. 


In our practical implementation, we used a Monte Carlo 
algorithm to create maps of the detectable rate of GRB 
per square degree per year from M31 and the Galactic fore¬ 
ground given a WFC pointing direction. The algorithm picks 
a burst in the Galaxy or M31 from the integrated density 
distribution between R — 0 and R — Ro and random an¬ 
gular coordinates. Then we compute Ry and the position 
on the sky, and from that the detection probability for the 
WFC. This probability is then added to the total in the ap¬ 
propriate sky location and/or flux bin and the procedure 
is repeated until the rate maps have sufficiently low Monte 
Carlo errors. In the case of beamed models, we also check 
whether the burst is shining in our direction before counting 
it as detectable. For this we make the approximation that 
GRB are ejected from a galaxy at such high velocities that 
they move at constant velocity in straight lines. We also as¬ 
sume they are all born at the centre of the galaxy and that 
their emitted flux is constant for all directions close enough 
to the direction of motion, and zero outside some critical an¬ 
gle. This simplification allows us to check the detectability 
of the burst using only the angle between our line of sight 
to the burst and the direction from the burst to the cen¬ 
tre of M31. (This approximation is quite good, because the 
neutron star formation rate drops exponentially from the 
centre with a small scale length of only about 4kpc, leading 
to an extra ‘smearing’ of the maps with an angular width 
of only Ssmear ~ 0.3°.) The opening angles (axis to edge) 
of the beaming cones we used for Rc =30, 40, 50, 60, and 
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Figure 5. Contours of the ex¬ 
pected number of bursts in 
a 2-month observation of un¬ 
beamed GRB (left) and a 6- 
month one of beamed GRB. 
The horizontal axis gives Rc, 
the vertical one the percent¬ 
age of GRB with an X- 
ray excess. The contour val¬ 
ues are 1 (thin-dotted), 2 (thin- 
dashed), 4 (thin-solid), 8 (thick- 
dotted), 16 (thick-dashed), and 
32 (thick-solid). Left is the 
number from M31, right that 
from the Galactic foreground. 


70 kpc are 20, 15, 12, 10, and 8 degrees, respectively, to keep 
RcO constant. Also, for these models we used outer radii of 
2 Mpc for the halos. (They should be unlimited, but beyond 
that virtually no objects will be detectable.) 


4 RESULTS 

There are three main variables on which the detectability of 
an excess of GRB to M31 depends most: the core radius, the 
fraction of bursts, /ex, with an X-ray excess, and beaming. 
To illustrate the dependence on Rc and beaming, we show 
in Fig.^ maps of the detectable rate of GRB as a function 
of position on the sky (in number per square degree per 
year). The top panels are for unbeamed bursts, the bottom 
ones for beamed bursts. The influence of Rc is reflected by 
the differences between the left (30 kpc) and right (60 kpc) 
panels. As noted above, the image becomes much smaller for 
beamed bursts and the size is independent of Rc, whereas 
for unbeamed bursts the increase of the image size with Rc 
is clear. Also note the increased peak rate (labelled ‘max’ in 
each panel) for smaller core radii. This is because the central 
density scales as R/®, and we do see a fair fraction of the 
bursts at the centre of M31. For comparison, the detectable 
rate from our own Milky Way is about 0.03/deg^/yr without 
beaming to 0.1/deg^/yr with beaming (the latter is greater 
due to the larger assumed outer radius for beamed models). 

Because beaming so much reduces the area of sky over 


which GRB can be seen, the total rate in the held of view 
will be much less for beamed bursts. This is very apparent 
in Fig. in which we show the expected number of bursts as 
a function of parameters. The top panels show the expected 
numbers in a 2-month WFC camera observation with M31 in 
the centre of the held of view for unbeamed bursts. Notice 
how the number from M31 always dominates strongly. In 
the beamed models (bottom panels), however, the numbers 
from both are comparable and a much longer observation is 
needed to detect an excess towards M31. 

Now we must create a practical test of whether seeing 
a certain number of GRB towards M31 constitutes evidence 
for or against a halo of GRB around it. Since the foreground 
bursts are more spread out on the sky than those of M31, 
this entails hnding the optimum area of the detector to use 
as the region within which we look for bursts. A good choice 
of boundary turns out to be a contour on which the summed 
rate of M31 and Milky Way bursts is a hxed fraction of the 
summed peak rate. For a given set of halo model parame¬ 
ters, we then hx an observing time and a boundary. This 
completely specihes the expected number of bursts, Emw, if 
there is no halo around M31, and the expected total number 
if there is, i?MW-i-M 3 i. The actual numbers we would get in 
an observation, Amw or Nmw-i-m31, have a Poisson distribu¬ 
tion around the expected values. Let us choose a threshold 
value Nth, which defines the boundary between the accept¬ 
ing the null hypothesis Ho ‘there is no halo with these pa- 
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Figure 

6. Contours of the probability 
of making the wrong decision 
based on a 1-month observation 
of unbeamed GRB (top) and a 
6-month one of beamed GRB. 
The horizontal axis gives Rc, 
the vertical one the percentage 
of GRB with an X-ray excess, 
/ex. The contour values range 
from 0.001 (thin-dashed) to 0.1 
(thick-solid) in steps of \/T0. 


rameters’ and accepting the alternative Hi ‘there is a halo 
with these parameters’. Obviously we can make two kinds 
of wrong decision: reject a halo when there is one, which 
happens with probability P{N < /fthli^MW-i-Msi) = Pi, and 
accept one if there is none, which happens with probability 
P{N > Althl^i'Mw) = P 2 - The best test will have small val¬ 
ues of Pi and P 2 ; Since we have no preference for either type 
of error above the other, we shall define A'th as the value for 
which Perr = max(Pi,P 2 ) is as small as possible. (The two 
can never be made equal because A'th can only assume inte¬ 
ger values.) For a given halo model and observing time we 
can then further optimise the test by varying the bound¬ 
ary of the detector region inside which we accept bursts to 
get the lowest overall Pen- In practice the optimal contour 
boundary lies at a rate of 5-10% of the peak rate (but the 
minimum is shallow). 

In Fig.|^ we show contours of optimised Perr values, 
again as a function of Rc and /ex. We can see that for the 
nominal value /ex = 0.5 we get a quite decisive test with 
Perr < 0.01 for all core radii if bursts are not beamed. The 
reason that the results do not depend as much on core ra¬ 
dius as one might have thought is that sampling distance is 
not the major issue: the X-ray brightest bursts can be seen 
out to a few Mpc, well beyond M31. An increase in core 
radius will decrease the luminosity of bursts, but it will also 
increase the central density of the halo and concentrate the 
bursts more near the detector centre. The net effect is not 
large. 

If GRB are beamed, Perr falls in the range 0.03-0.1 even 
for a 6-month observation, and the case is rather less con¬ 
vincing. Should the bursts be strongly beamed, our results 
point to a very cheap (in space dollars) and useful satellite 
mission that can settle the issue: a camera similar to the 
WFC, but with half the held of view and therefore only one 
quarter of the background. It should be pointed to M31 for 
3-6 months and then in some other direction for a similar 
period of time and would either detect a corona of beamed 
sources or rule it out. Since many known X-ray sources in 
M31 would be easily detected in a small fraction of the re¬ 
quired observing time by the same instrument, there would 
be considerable beneht to such a mission for the study of 


variability and population analysis of bright X-ray binaries 
in M31. 


5 DISCUSSION AND CONCLUSION 

Previous work on constraining coronal models of GRB has 
many similarities with our own calculations. Liang (1991) 
found that ROSAT might detect some GRB in X rays as¬ 
suming a now abandoned disc model for the distribution of 
GRB. Li, Fenimore, & Liang (1996) used a similar method 
to our own both for beamed and unbeamed models. Their 
calculations differ from ours mostly in that they use neither 
the spread in X-ray luminosity nor the X-ray excesses. Their 
hypothetical instrument had Smin = 0.1 phot cm“^ s“^. As 
we can see from Fig.^ 30-90% of the bursts that the WFG 
can see on-axis are below this limit, so it is no surprise that 
we find more optimistic prospects for detecting GRB in M31. 

Harrison and Thorsett (1996) considered a variety of 
real instruments, calculating the detectable rate in much 
the same way as we did (including the spectral variability 
using the same set of spectra from Band et al. 1993). They 
conclude that only a novel instrument sensitive to photons 
in the 10-200 keV range and with a field of view of 18° would 
be capable of detecting M31 in one year. While they did not 
include the possibility of X-ray excesses, they would without 
doubt have realised the potential of the SAX WFG if they 
had included them in their work. 

In summary, we have shown that the hitherto neglected 
spread in X-ray to gamma-ray luminosity ratios of gamma- 
ray bursts substantially increases the prospects for deciding 
the gamma-ray burst distance scale. The case is further im¬ 
proved greatly by the recent discovery that a substantial 
fraction of gamma-ray bursts have X-ray excesses (Preece 
et al. 1996). A one-month observation of M31 with an ex¬ 
isting instrument, the SAX WFG, will be decisive for es¬ 
tablishing whether or not the Andromeda Nebula harbours 
a population of bursters, unless bursters only emit radia¬ 
tion in fairly narrow cones along their direction of motion. 
In that case, a dedicated, cheap mission similar to the SAX 
WFC should resolve the issue in about one year of observing 
time. Observing proposals to do the experiment in WFC sec- 
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ondary (i.e. unguaranteed) time have been accepted, so the 
gamma-ray burst distance scale may not remain uncertain 
much longer. 
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